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We study an exciton in a type II quantum dot, where 
the electron is confined in the dot, but the hole is located in 
the barrier material. The exciton properties are studied as 
a function of a perpendicular magnetic field using a Hartree- 
Fock mesh calculation. Our model system consists of a pla- 
nar quantum disk. Angular momentum (Z) transitions are 
predicted with increasing magnetic field. We also study the 
transition from a type I to a type II quantum dot which is 
induced by changing the confinement potential of the hole. 
For sufficiently large magnetic fields a re-entrant behaviour is 
found from lh = to lh 7^ and back to lh = 0, which results 
in a transition from type II to type I. 

PACS: 73.21.La, 71.35.Ji, 85.35.Be 



I. INTRODUCTION 

Self-assembled quantum dots Jl]] have become the sub- 
ject of intensive research, both theoretically and experi- 
mentally, since their first realization in the early nineties 
The reason for this large interest is e.g. due to their 
possible applications in opto-electronic devices, such as 
quantum dot lasers. The formation of this type of dots by 
the Stranski-Krastanow growth mode requires two semi- 
conductor materials with a considerable lattice mismatch 
of typically 5%. Many experimental ]6| |lO|| and theoret- 
ical [Il|-[l5[ works are devoted to type I structures, e.g. 
InAs/GaAs or InAlAs/AlGaAs, where both electrons 
and holes are located inside the quantum dots. 

Also very interesting, though yet less studied, are the 
type II quantum dots, where the quantum dot forms an 
antidot for one of the types of carriers, e.g. for the holes 
in typically the InP/GalnP system or the electrons in 
e.g. GaSb/GaAs. Landau level formation in strongly 
optically populated type II dots was observed by No- 
mura et al. |l6| in the photoluminescence spectra at high 
magnetic fields. Other magneto-photoluminescence ex- 
periments on vertically stacked InP quantum dots were 
performed by Hayne et al. p7| . Sugisaki et al. M] stud- 
ied the magnetic field effects in a single InP dot. The 
optical recombination spectrum and the carrier dynamics 
of the GaSb/GaAs system have been studied experimen- 
tally by Hatami et al. |L9| . 

Whereas the type I system has been the subject of 
many theoretical works, only few theoretical studies have 
paid attention to the type II system. Pryor et al. pQ| 
studied the electronic structure of InP/GalnP, using a 
strain-dependent kp Hamiltonian. Also for InP/GalnP 



dots, Nomura et al. (2TJ performed a theoretical calcu- 
lation of the Landau levels in a high magnetic field, by 
solving the Hartree equations self-consistently. Using the 
Hartree-Fock approximation, the binding energy of exci- 
tons, charged excitons and biexcitons was studied by Le- 
long et al. |E| in GaSb/GaAs dots at zero magnetic field. 
The magneto-exciton in a GaSb/GaAs dot was investi- 
gated by Kalameitsev et al. |23| . They found transitions 
of the angular momentum with increasing magnetic field. 

In the present paper, we focus our attention to the 
properties of a single exciton, which is bound by the 
Coulomb interaction in a model type II quantum dot. We 
will take material parameters of the InP/ GalnP system. 
Furthermore, we apply an external magnetic field in the 
growth direction, i.e. B = Be z . Including a magnetic 
field allows us to investigate the transition region from 
exciton confinement due to the Coulomb potential, to a 
confinement which is due to the magnetic field. As a 
model system we take a planar quantum disk, and as- 
sume that the particles are confined in a plane in the 
z-direction. Strain effects are neglected in this model 
system. 

In our model type II quantum dot the electron is con- 
fined in the dot and the hole sits outside. The corre- 
sponding geometry is shown in Fig. 1. The reverse con- 
finement situation will lead to the same physics. As we 
do not take the confinement effects due to strain into ac- 
count, the hole is only confined because of the Coulomb 
attraction to the electron. As we have no a priori knowl- 
edge about the width of the hole wavefunction it is dif- 
ficult to choose good basisfunctions for the expansion of 
the hole wavefunction. Therefore, we solved the Hartree- 
Fock (HF) equations on a grid, which allows very flexible 
solutions, in principle of arbitrary shape. With the same 
motivation, similar Hartree-Fock mesh calculations were 
recently used in atomic physics p3]. 

As confinement potential, we take hard walls of finite 
height. By varying the hole confinement potential, we 
can study the transition from a type I structure (i.e., the 
hole is confined in the dot) to a type II structure (i.e., the 
dot is a barrier for the hole). We show that for small anti- 
dots the attraction of the hole to the electron is stronger 
than the barrier energy and the system is still type I. 
Increasing the barrier height and/or the size of the dot 
induces a transition to a type II system. Furthermore, 
we found angular momentum transitions with increas- 
ing magnetic field. For large enough magnetic fields (de- 
pending on the height of the potential barrier) , we find a 
new re-entrant behaviour to the zero angular momentum 
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state. 

The paper is organized as follows. In Sec. II, we de- 
scribe briefly our theoretical model. The numerical re- 
sults are presented in Sec. III. In Part A of this Sec- 
tion, we discuss the effect of a varying magnetic field and 
explain the origin of the angular momentum transitions. 
Part B deals with the transition from a type I to a type II 
system. Part C is dedicated to the re-entrant behaviour. 
In the last Part, D, we present the results for the excita- 
tion spectrum. Our results are summarized in Sec. IV. 
In the Appendix, we discuss in more detail the method 
we used for the calculation of the Hartree integral. 



II. THEORETICAL MODEL 
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where K(x) is the complete elliptic integral of the first 
kind. More details about the calculation and numerical 
implementation of this integral is given in the Appendix. 
After convergence, the total energy is given by 
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The contribution of the correlation to the total energy 
is neglected in HF, but for the self-assembled quantum 
dots, it is expected to be less than 2% [Q. 



The energies and wavefunctions are obtained by solv- 
ing the following HF single particle equations in the ef- 
fective mass approximation (with m e and rrih the ef- 
fective electron and hole masses, respectively, r e h = 

= eB/m e and uj c ,h — eB/rrih) 
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III. RESULTS 

A. Angular momentum transitions 

First, we calculated the groundstate energy of the exci- 
ton as a function of the external magnetic field. We took 
the following parameters: m e = 0.077to , — 0.6to , 
V e = 250meV, Vh — —50meV, and e = 12.61, which are 
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where we made use of the axial symmetry by taking 
f e (r e , <p e ) = e ll ^^ e (r e ) and V h (r h , <p h ) = e lt ^^ h (r h ) 7 
and where the densities p e (f') and Phir') are given 
by respectively \^ e (r e , Lp e )\ 2 and \^ h {r h ,Lp h )\ 2 . The 
Hartree-fock equations were solved using a finite differ- 
ence scheme. More details about the implementation of 
this finite difference scheme can be found in Refs. |^5|,|l5| . 
Note that there are no exchange terms as we only con- 
sider a single electron and a single hole. However, these 
equations can still be called HF as the self-interaction is 
excluded. As confinement potentials we take hard walls 
of finite height: 



V e ,h(r e ,r h ) 



0, otherwise, 



(2) 



with R the radius of the disk, and where we took V e pos- 
itive and Vh negative. Note that the only good quantum 
number is the total angular momentum in the z-direction, 
defined by L = l e + 1^. 

These equations must be solved self-consistently, which 
is done iteratively. We start with the free electron solu- 
tion because in the absence of any Coulomb interaction 
only the free electron is confined. The Hartree integrals 
are integrated numerically 



typical for the InP/GalnP system |20 and consider a 
dot of radius R = 8nm |17[| . Our numerical results are 
d^pftked in Fig. 2 and show that the exciton groundstate 
exhibits transitions of the angular momentum lh as a 
function of the magnetic field (indicated by the arrows). 
These changes in angular momentum of the groundstate 
are not present in type I dots and are a direct conse- 
q{iifeibe of the fact that we are dealing with type II dots. 
For the hole, the disk acts as a barrier and by increasing 
the magnetic field the hole is pushed closer to the disk 
boundary which leads to an increase of the hole potential 
energy. For a certain magnetic field it is for the hole en- 
ergetically more favourable to jump to a higher lh state, 
which brings the hole further away from the disk inter- 
face. This is also demonstrated in the inset of Fig. 2, 
where a contourplot of the density of the hole wavefunc- 
tion is shown, as a function of both the magnetic field 
and the radial position. It is apparent that the hole is 
located close to the disk, even at zero magnetic field and 
that with increasing magnetic field the hole is pushed 
closer to the border of the disk. At the angular momen- 
tum transitions the hole is spread out a little more, and 
jumps a distance away from the disk interface. But note 
that on the average the hole is pushed closer to the disk 
boundary and its width decreases with increasing mag- 
netic field. For the present case we find five transitions 
for a magnetic field up to B = 50T. 

The magnetic field values at which the angular mo- 
mentum transitions occur will depend on the disk radius 
R. Fig. 3 shows a phase diagram of the lh transitions as 
a function of the magnetic field B and the disk radius 
R, for the InP parameters used above. With increasing 
disk radius, the transitions shift to lower magnetic field 
values. This can be understood as follows: for larger 
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disks, a smaller magnetic field is needed to push the hole 
close to the border of the disk, and to induce an angular 
momentum transition. 

Another way to understand the angular momentum 
transitions is as follows. The hole is spatially confined 
into a ring-like area, and if we make the extreme sim- 
plification of a zero width ring, the hole energy is given 
by 



2mR 2 



(5) 



from which it is clear that the ground state exhibits angu- 
lar momentum transitions each time the flux through the 
ring <j) equals (lh + l/2)0o with 0o = hc/e the quantum 
of flux. 



B. Type II to type I transitions 

In Fig. 4 we give a closer look at the hole wavefunc- 
tion for a very small disk R — 2nm, Vh — — 50meV and 
B = OT, and we find that, even without a magnetic field, 
the hole is partially situated inside the quantum disk. 
This is a remarkable effect, as from the shape of the ef- 
fective potential, defined by the sum of the confinement 
potential and the Hartree potential (inset of Fig. 4), we 
would expect the hole wavefunction to be situated in the 
barrier. We attribute this effect to a kind of tunneling of 
the hole through the quantum disk, as a consequence of 
the very small disk radius. Remark that this state con- 
tinues to be the groundstate and that the effect becomes 
even stronger for higher fields. Another consequence of 
this effect is a higher overlap of the electron and hole 
wavefunctions, which is an indication of type I behaviour. 

In a next step, we studied the exciton properties as a 
function of the hole confinement potential, which allows 
us to explore the transition region from type I systems 
(Vh > 0) to type II systems (Vh < 0) . Hereby we kept the 
disk radius fixed at R = 8nm. Fig. 5 shows the phase dia- 
gram for the angular momentum transitions as a function 
of the confinement potential Vh and the magnetic field B. 
A feature that immediately catches the eye is that up to 
Vh — —24.5meV the lh = state remains the groundstate 
over the total i?-region under consideration. Investigat- 
ing this more in depth, we find that, even for B = OT, 
the hole wavefunction is located almost entirely inside 
the quantum disk. This is a consequence of the Hartree 
potential (due to the attraction to the electron) which 
overcomes the potential barrier of the disk. Therefore, 
we can speak of type I systems up to Vh — ~24.5meV. 

In order to have a more physical idea of the origin the 
type I to type II transition, we developed the following in- 
tuitive picture. When the system is type I, the hole will 
be located inside the quantum disk. This will happen 
when the effective potential (the sum of the hole confine- 
ment potential and the Hartree potential, i.e. V c (rh) + Vh 
(see the inset of Fig. 4)) is lower at = than at r/, = R. 



For a type II system however, the effective potential will 
be lower at the boundary, and the hole will prefer to sit 
outside the disk but near the boundary. Generally, we 
can state that the transition from type I to type II oc- 
curs when the effective potential at the origin equals the 
one at the radial boundary. Because we take the con- 
finement potential zero inside the disk, this leads to the 
following formula: 



V c (r h = 0) = V h + V c (r h = R). 



(6) 



From this equality, we can find an estimate for the con- 
finement potential Vh at which the type I to type II 
transition occurs. The Hartree potential was calculated 
within the approximation of an infinitely high electron 
confinement potential, where the electron wavefunction 
can be expressed by the Bessel functions Jo, i.e. 
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with N 2 = 2/(R 2 \J 1 (kr' e )\ 2 ) the normalization, and 
k = \/2m e E/h. In Fig. 6, we plot the values for Vh 
as a function of the disk radius R at which such a type 
I to type II transition occurs. We show the results ob- 
tained by the full Hartree-Fock calculation (full curve, 
squared dots) and the approximated results, as obtained 
from Eq. (6) (dashed curve, circular dots). We find that 
for large R the two curves converge to each other, but 
that for small R an appreciable discrepancy exists. This 
is due to the fact that for small disk radii, the approx- 
imation of a hard wall confinement is less justified. By 
pushing the electron wavefunction completely into the 
disk (in contrast to the 'real' case of a finite potential, 
where the wavefunction can tunnel into the barrier), the 
Hartree potential is strongly enhanced, thereby leading 
to a strong enhancement of the critical confinement po- 
tential Vh- 

Another interesting property is the probability for re- 
combination of the exciton. This is proportional to the 
square of the overlap integral 23 
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Notice that the first integral is equal to 2TrSi ii+ i h , what 
means that the probability for de-excitation is only non- 
zero for the case l e + lh = 0. This implies that after an 
angular momentum transition the probability for recom- 
bination of an exciton decreases drastically. In photolu- 
minescence (PL) experiments, one will observe a strong 
quenching or even disappearance of the PL spectrum af- 
ter a certain value of the magnetic field. 

Fig. 7 shows the overlap integral I as a function of 
Vh, for B = OT. Without a magnetic field, the lh = 
state is always the groundstate and therefore, / is non- 
zero over the total region. Up to Vh = ~25meV the 
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overlap is large, and further increasing — Vh, we find a 
sudden strong decrease of the overlap. The reason for 
this behaviour is directly related to the position of the 
hole wavefunction. As long as the hole is sitting inside 
the disk, the overlap will be very large. However, from 
the moment the hole jumps outside the disk, the over- 
lap decreases strongly. In fact, from Fig. 7 we can infer 
immediately the position of the hole. Furthermore, the 
region of the strong decrease in overlap indicates the tran- 
sition from type I to type II behaviour. The dashed line 
gives the overlap integral at B = 2QT. As we see from 
Fig. 5, a transition to the 1^ = 1 state occurs when the 
confinement potential Vh approaches —2A.5meV and the 
condition l e + lh = for recombination of the exciton, is 
no longer satisfied, which leads to I = 0. The recombi- 
nation of the exciton will happen through indirect pro- 
cesses, resulting in a much longer lifetime of the exciton. 
The changing lifetime can be detected experimentally by 
a changing lineshape |27]j2q] . 



C. Re-entrant behaviour 

In this section, we concentrate more closely on the type 
I - type II transition region, i.e. Vh between —20meV 
and — 30meV. As an example, we investigated the exci- 
ton groundstate energy for Vh = —27meV. The result is 
depicted in Fig. 8 which shows one additional remarkable 
feature: after several lh transitions with increasing mag- 
netic field, we find at sufficiently large magnetic field, i.e. 
B ~ 42T, a re-entrance of the lh = state. It is interest- 
ing also to take a look at the evolution of the wavefunc- 
tion with increasing magnetic field. This is depicted as a 
contourplot in the inset of Fig. 8. Initially, for very small 
magnetic fields, we find that a small part of the wave- 
function has already entered the dot region. However, at 
B ~ 6T, due to a jump to a higher angular momentum 
state, the hole wavefunction is pushed outside the dot re- 
gion. Further increasing the magnetic field leads to more 
lh transitions, as we found already in Section III. A. At 
the specific magnetic field value however, where the re- 
entrance of lh = occurs, we find that suddenly the hole 
wavefunction jumps almost entirely inside the disk. At 
this point, the magnetic field and the attraction of the 
electron overcomes the potential barrier of the quantum 
disk, and it will be energetically more favourable for the 
hole to sit inside the disk. 

The re-entrant behaviour is also visible in the [B, Vhj- 
phase diagram (Fig. 5). We want to emphasize that 
there will be a re-entrant behaviour for any value of 
Vh < —24.5meV, for sufficiently large magnetic fields. 
This can already be seen from Fig. 5, where the line 
which indicates the transition from a certain Z^-state to 
the lh = state is not a straight vertical line, but has 
a small slope. For example for Vh = — 50meV we found 
that a magnetic field of B = 193T is needed to induce 
this re-entrance to the lh = state. 



Another question which arises is how the disk radius in- 
fluences the re-entrant behaviour. In our previous inves- 
tigation of the influence of the disk radius on the ground- 
state energy, we found no evidence of this, because we 
did not consider large enough magnetic fields for the con- 
finement potential under consideration (Vh = —50meV) . 
Therefore, we decided to make a new (B, i?)-phase dia- 
gram, this time for Vh = — 25meV, for which we know 
from Fig. 5 that re-entrant behaviour occurs at rather 
small magnetic fields. The result is depicted in Fig. 9, 
and the first striking feature is that the re-entrant be- 
haviour occurs for any disk with radius R > 8nm. Fur- 
thermore, we find more lh transitions for larger disk radii, 
hereby increasing the magnetic field at which the re- 
entrance of lh = takes place. However at R ~ 12nm, we 
find that the magnetic field position of the re-entrant be- 
haviour reaches a maximum value, B ~ 42T. For larger 
disk radii, the re-entrance occurs at slightly decreasing 
magnetic field. Indeed, with increasing disk radius, elec- 
tron and hole are drawn more and more apart, and there- 
fore it will sooner become energetically more favourable 
for the hole to jump inside the disk. 

This re-entrant behaviour can be understood qualita- 
tively from the following simple model. We compare the 
approximate energies for a hole of respectively a type II 
and a type I system: 
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The approximate energy for a hole in a type II system 
(Eq. (9a)) is constructed by approximating the disk (with 
radius R) by an infinitely high barrier for the hole. The 
first term of Eq. (9a) gives the one-particle energy of 
the hole, which is just the first energy level of a particle 
in a magnetic field which fulfils the zero-wavefunction 
condition at the disk boundary. The second term is the 
Coulomb interaction between the electron and the hole, 
and the third is the potential energy of the hole. In a type 
I system (Eq. (9b)) both electron and hole are located 
inside the disk, and are subjected to a magnetic field. 
Now the one-particle hole energy is just the energy of the 
first Landau level. We approximated the Coulomb energy 
(second term in Eq. (9b)) by using the single particle 
wavefunction of electron and hole. At strong magnetic 
fields and for large disk radii i.e. R > Ib, the potential 
confinement by the disk can be neglected with regard to 
the confinement by the magnetic field. The single particle 
wavefunctions are then the well-known wavefunctions of 
a particle in a magnetic field, given by 
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for the groundstate (n = 0, lh = 0). As the wavefunctions 
are mass- independent, there is no distinction between the 
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electron and the hole. Therefore we can treat our system 
as being completely analogous to a system consisting of 
two electrons in a magnetic field, for which our first order 
approximation of the Coulomb interaction energy can be 
calculated analytically, and which reduces to 
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Equalization of Eqs. (9a) and (9b) gives us the mag- 
netic field at which the re-entrance occurs as a function 
of both disk radius R and confinement potential Vh ■ The 
transition magnetic field can be obtained analytically as 
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For a fixed radius R, we can vary Vh and deduce the 
magnetic field at which the transition occurs. The result 
is shown by the dashed curve in Fig. 5. For large fields, 
the approximated curve (dashed line) has qualitatively 
the same behaviour as the curve obtained from the full 
Hartree-Fock treatment and the two curves converge to 
each other. The result for a fixed Vh, when varying the 
disk radius R, is shown by the dashed curve in Fig 9. 
We find a perfect agreement for large disk radii, where 
our model is valid. The discrepancy for small disk radii 
is a consequence of our assumption to neglect the disk 
confinement. 

Fig. 10 shows the overlap integral J as a function of 
the magnetic field, for confinement potentials of the hole 
Vh = —25meV (solid curve) and Vh = —27meV (dashed 
curve) and a fixed disk radius R = 8nm. At first we find 
a slowly increasing overlap, which is a consequence of the 
increasing magnetic field, pushing the particles closer to- 
gether. The already rather large value of the overlap in- 
dicates that a considerable part of the hole is already 
situated inside the disk. When the first /^-transition 
occurs, the overlap falls immediately down to zero, be- 
cause l e + lh 7^ 0. The overlap remains zero, until the 
lh = state returns as the groundstate and the condi- 
tion l e + lh — is satisfied. This re-entrance of lh = is 
accompanied by a jump of the wavefunction in the disk, 
and this leads to the strong enhancement of the overlap 
value. We see that the re-entrance of the lh = state 
happens at lower magnetic fields for the lower potential 
barrier. 

Fig. 11 shows the overlap integral I as a function of 
the disk radius R for B = 0T (solid curve), B = 20T 
(dashed curve) and B = 40T (dotted curve). This figure 
gives evidence for the fact that for very small radii the 
hole wavefunction is almost entirely situated inside the 
disk. For increasing disk radius, the hole is pushed more 
and more outside the disk, hereby decreasing the value 



of the overlap integral. This decrease is initially less for 
increasing magnetic field because of the enhanced local- 
ization effect. For sufficiently large magnetic fields, lh 
transitions are induced, which leads to a zero overlap 
integral. Also here we see that for sufficient large R a re- 
entrant behaviour to the lh = state is found at which 
point the overlap integral becomes again non-zero. 



D. Excitation spectrum 

Lastly, we investigated the exciton energy spectrum as 
a function of the magnetic field. The physical parame- 
ters used in this calculation are the ones mentioned in 
Section II, with a disk radius R — 8nm. We considered 
states with different radial quantum numbers k e , kh and 
ifferent angular momenta l e , lh- Note that k e and kh are 
pproximate quantum numbers for the current system. 
^Furthermore, our results are a first order perturbation 
theory approximation to the real energy spectrum, as we 
perform only one HF iteration, i.e. the energy of the ex- 
citon is obtained by solving the equation for the hole in 
the field of the confined electron. Note that changing the 
electron quantum number results in a strong increase of 
the energy value. The states (k e ,l e ) = (2,0) and (1,±2) 
appear to be already unbound, i.e. the energy exceeds 
the electron barrier of 250meV. The inset of Fig. 12(a) 
shows the bound states of the energy spectrum where we 
varied both k e and l e , keeping kh and lh fixed at (0,0) . 
The main part of Fig. 12(a) shows the energy spectrum 
for fixed (k e ,l e ) = (0,0) and varying the hole quantum 
numbers kh and lh- We find that now the energy values 
span a smaller energy region. This is due to the fact 
that: i) the hole is much heavier than the electron, and 
therefore has substantially lower energies, and ii) the hole 
is less confined. In fact, for every possible value of the 
electron quantum numbers k e and l e , one has a spectrum 
of all possible (kh,lh) values, and because these span a 
smaller energy region, the total energy spectrum will con- 
sist of mainly the electron branches, with superimposed 
on each of them the spectrum with the changing hole 
quantum numbers. 

Notice also the anomalous behaviour of certain states, 
e.g. (khjh) = (2,0) in the high magnetic field region. 
To investigate this further, we concentrated on the vari- 
ation of the radial quantum number kh, keeping the an- 
gular momentum lh fixed at 0. This result is shown in 
Fig. 12(b) which clearly shows the occurrence of anti- 
crossings. These anti-crossings are due to the fact that 
the radial quantum number is not a good quantum num- 
ber, leading to strong mixing of radial states at the anti- 
crossings. 
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IV. CONCLUSIONS 

We investigated the exciton properties in a strongly 
simplified type II model quantum disk, with the hole lo- 
cated in the barrier. Strain effects were disregarded and 
a flat disk geometry was assumed. Because in our model 
system there is no geometrical confinement for the hole, 
the only "confinement" comes from the attraction to the 
electron, i.e. the Coulomb interaction energy. We solved 
this problem by using a Hartree-Fock mesh calculation, 
which allowed us to calculate the exciton energy, without 
an a priori knowledge of the single particle hole wavefunc- 
tion. 

We studied the influence of a perpendicular applied 
magnetic field, and found angular momentum transitions 
with increasing magnetic field. These are a consequence 
of the fact that the magnetic field pushes the hole closer 
to the disk, making it energetically more favourable to 
jump to a higher lh state. Varying the disk radius showed 
that the transitions shift to lower magnetic field for larger 
R. We also found that the hole is located almost entirely 
inside the disk for very small disk radii. 

Furthermore we investigated the transition region from 
type I to type II systems, by varying the confinement 
potential of the hole, Vh- A striking feature here is the 
fact that we are dealing with type I systems up to Vh > 
— 24.5meV. Even at B = QT, the Coulomb attraction 
overcomes the potential barrier and the hole is situated 
inside the quantum disk. Taking a closer look at the 
transition region between type I and type II, showed the 
existence of a re-entrant behaviour of the lh = state. 
This re-entrant behaviour is coupled with a sudden jump 
of the wavefunction into the disk. 

The angular momentum transitions and the re-entrant 
behaviour should be measurable experimentally by 
quenching of luminescence and/or changing lineshapes. 

In a last part, we studied the excitation spectrum as 
a function of the magnetic field. We varied the quantum 
numbers k and I for both electron and hole, and found 
that for every value of (k e ,l e ) one has a spectrum con- 
sisting of the different radial and angular momentum hole 
states. Furthermore, taking a closer look to the varying 
kh states, with fixed lh and fixed electron quantum num- 
bers, we found an anti-crossing of levels, a consequence 
of the fact that kh is not a good quantum number and 
therefore lifts the degeneracy. 
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VI. APPENDIX: CALCULATION OF THE 
HARTREE INTEGRAL 

The Hartree integral, which expresses the effect of one 
particle, e.g. the electron, on the other particle, e.g. the 
hole, is given by 



J - J *><> J « 



„(/) 



yjr 2 + r' 2 — 2rr' cos(ip — ip') 
(Al) 



As we arc dealing with cylindrical symmetry, we can 
remove the (^-dependence. The integral over the angle 
becomes the complete elliptic integral of the first kind, 
which converts (Al) into 



p(r')r' ( 4rr' 



4 / r -^—K ( . " ' ,.. ) dr'. 
r + r' \(r + r') 2 / 



(A2) 



The radial integral has to be solved numerically. We use 
a polynomial approximation for the elliptic function [fjjif , 
namely 

K(x) = [a + aix' + a 2 x' 2 } - [b + bix' + b 2 x' 2 ] ln(x'), 

(A3) 

with x' = 1 — x and where the coefficients at and bt are 
given in Ref. p(|. Since this implies the appearance of 
a logaritmic divergence in the integrand, the commonly 
used trapezoidal rule will give bad results. Therefore 
we used another method, the so-called 'logaritmically 
weighted method' which takes into account this problem. 
Generally, the following integral can be considered: 



I(r) = l 1 d xF(x)ln^ > 
which, after transformation, becomes 



(A4) 



iV-l -h 

I(r) = V / dxF(x + hi) In 



{x+{r + hi)) 2 ' 1 ' 



with h the discretization step and N the number of steps. 
If we replace F(x + hi) by + [Fi+i — Fi\ (x/h) , we can 
write (A5) as 



JV-l 



I(r) = £ {FMr) + [F i+1 - Fildir)}, (A6) 



i=0 



and the remaining problem is the calculation of the coef- 
ficients Ai(r) and Ci(r). The integrals which determine 



G 



the coefficients can be solved exactly, which leads to the 
following results: 

Ai(r) = a(h - (r - hi)) - a{h + (r + hi)), (A7) 

with 

a(y) = 2y\ny + 2(h-y)ln(h-y), (A8) 

and 

C\(r) = h-^cih - (r - hi)) - c(h + (r + hi))} - 2r, 

(A9) 

with 

c(y) = (y(2h - y)) lny + (h- yf In (h - y) . (A10) 
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FIG. 1. The geometry of the system under consideration, 
with side and top view. 

FIG. 2. The exciton energy as a function of the magnetic 
field, for R = 8nm, V e = 250meU and V h = -50meV. The 
successive /^-transitions are indicated by arrows. The inset 
shows a contourplot of the square of the hole wavefunction as 
a function of the magnetic field. 

FIG. 3. Phase diagram of the In transitions as a function 
of magnetic field B and disk radius R, for Vh = — 50meU. 

FIG. 4. Hole wavefunction for R — 2nm and 
Vh — —50meV, at B = 0T. The inset shows the effective 
confinement potential for the hole. 
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FIG. 5. Phase diagram of the successive l h states for vary- 
ing confinement potential 14, as a function of the magnetic 
field. The disk radius is fixed at R — 8nm. The dashed line 
indicates the result obtained by the approximate model. 

FIG. 6. Confinement potential of the hole at which a tran- 
sition from a type I to a type II system occurs, as a function of 
the disk radius R. The solid curve (squared dots) indicates the 
result obtained within the HF treatment, whereas the dashed 
curve (circular dots) the approximated result is. 

FIG. 7. Overlap integral for a varying confinement poten- 
tial of the hole, for R — 8nm and at B — OT (solid curve) 
and B = 20T (dashed curve). 

FIG. 8. Exciton groundstate energy for Vh = —27meV. 
The successive lh states are indicated by arrows. Notice the 
re-entrant behaviour of the lh = state. 

FIG. 9. (B,R) Phasediagram for V h = -25meV. The 
re-entrant behaviour slowly decreases for large R. The dashed 
curve indicates the result obtained by the approximate model. 

FIG. 10. Overlap integral as a function of the magnetic 
field, for fixed R = 8nm and for Vh = — 25meV (solid curve) 
and Vh = —27meV (dashed curve). When l e + h =fc 0, the 
overlap is 0. 

FIG. 11. Overlap integral as a function of the disk radius, 
for Vh = — 25meV and magnetic fields of respectively OT 
(solid curve), 20T (dashed curve) and 40T (dotted curve). 



FIG. 12. (a) Energy spectrum for different kh and lh, 
for fixed (k e ,l e ) = (0,0), as a function of the magnetic 
field. Inset: idem, but now for different k e and l e , and for 
(kh,lh) = (0,0). (b) The same as in (a), but now for lh = 0. 
Notice the anti-crossings as a function of the magnetic field, 
which is due to the lifting of the degeneracy, or a strong mix- 
ing of the radial states. 
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